Spectral plots and the representation and interpretation of biological data

It is basic question in biology and other fields to identify the char- acteristic properties that on one hand are shared by structures from a particular realm, like gene regulation, protein-protein interaction or neu- ral networks or foodwebs, and that on the other hand distinguish them from other structures. We introduce and apply a general method, based on the spectrum of the normalized graph Laplacian, that yields repre- sentations, the spectral plots, that allow us to find and visualize such properties systematically. We present such visualizations for a wide range of biological networks and compare them with those for networks derived from theoretical schemes. The differences that we find are quite striking and suggest that the search for universal properties of biological networks should be complemented by an understanding of more specific features of biological organization principles at different scales.Comment: 15 pages, 7 figure

Kernel reconstruction for delayed neural field equations

Understanding the neural field activity for realistic living systems is a challenging task in contemporary neuroscience. Neural fields have been studied and developed theoretically and numerically with considerable success over the past four decades. However, to make effective use of such models, we need to identify their constituents in practical systems. This includes the determination of model parameters and in particular the reconstruction of the underlying effective connectivity in biological tissues. In this work, we provide an integral equation approach to the reconstruction of the neural connectivity in the case where the neural activity is governed by a delay neural field equation. As preparation, we study the solution of the direct problem based on the Banach fixed point theorem. Then we reformulate the inverse problem into a family of integral equations of the first kind. This equation will be vector valued when several neural activity trajectories are taken as input for the inverse problem. We employ spectral regularization techniques for its stable solution. A sensitivity analysis of the regularized kernel reconstruction with respect to the input signal u is carried out, investigating the Frechet differentiability of the kernel with respect to the signal. Finally, we use numerical examples to show the feasibility of the approach for kernel reconstruction, including numerical sensitivity tests, which show that the integral equation approach is a very stable and promising approach for practical computational neuroscience

Modeling Brain Resonance Phenomena Using a Neural Mass Model

Stimulation with rhythmic light flicker (photic driving) plays an important role in the diagnosis of schizophrenia, mood disorder, migraine, and epilepsy. In particular, the adjustment of spontaneous brain rhythms to the stimulus frequency (entrainment) is used to assess the functional flexibility of the brain. We aim to gain deeper understanding of the mechanisms underlying this technique and to predict the effects of stimulus frequency and intensity. For this purpose, a modified Jansen and Rit neural mass model (NMM) of a cortical circuit is used. This mean field model has been designed to strike a balance between mathematical simplicity and biological plausibility. We reproduced the entrainment phenomenon observed in EEG during a photic driving experiment. More generally, we demonstrate that such a single area model can already yield very complex dynamics, including chaos, for biologically plausible parameter ranges. We chart the entire parameter space by means of characteristic Lyapunov spectra and Kaplan-Yorke dimension as well as time series and power spectra. Rhythmic and chaotic brain states were found virtually next to each other, such that small parameter changes can give rise to switching from one to another. Strikingly, this characteristic pattern of unpredictability generated by the model was matched to the experimental data with reasonable accuracy. These findings confirm that the NMM is a useful model of brain dynamics during photic driving. In this context, it can be used to study the mechanisms of, for example, perception and epileptic seizure generation. In particular, it enabled us to make predictions regarding the stimulus amplitude in further experiments for improving the entrainment effect

Balancing with Vibration: A Prelude for “Drift and Act” Balance Control

Stick balancing at the fingertip is a powerful paradigm for the study of the control of human balance. Here we show that the mean stick balancing time is increased by about two-fold when a subject stands on a vibrating platform that produces vertical vibrations at the fingertip (0.001 m, 15–50 Hz). High speed motion capture measurements in three dimensions demonstrate that vibration does not shorten the neural latency for stick balancing or change the distribution of the changes in speed made by the fingertip during stick balancing, but does decrease the amplitude of the fluctuations in the relative positions of the fingertip and the tip of the stick in the horizontal plane, A(x,y). The findings are interpreted in terms of a time-delayed “drift and act” control mechanism in which controlling movements are made only when controlled variables exceed a threshold, i.e. the stick survival time measures the time to cross a threshold. The amplitude of the oscillations produced by this mechanism can be decreased by parametric excitation. It is shown that a plot of the logarithm of the vibration-induced increase in stick balancing skill, a measure of the mean first passage time, versus the standard deviation of the A(x,y) fluctuations, a measure of the distance to the threshold, is linear as expected for the times to cross a threshold in a stochastic dynamical system. These observations suggest that the balanced state represents a complex time–dependent state which is situated in a basin of attraction that is of the same order of size. The fact that vibration amplitude can benefit balance control raises the possibility of minimizing risk of falling through appropriate changes in the design of footwear and roughness of the walking surfaces

Search for the pair production of light top squarks in the e(+/-)mu(-/+) final state in proton-proton collisions at root s=13 TeV

A search for the production of a pair of top squarks at the LHC is presented. This search targets a region of parameter space where the kinematics of top squark pair production and top quark pair production are very similar, because of the mass difference between the top squark and the neutralino being close to the top quark mass. The search is performed with 35.9 fb(-1) of proton-proton collisions at a centre-of-mass energy of root s = 13 TeV, collected by the CMS detector in 2016, using events containing one electron-muon pair with opposite charge. The search is based on a precise estimate of the top quark pair background, and the use of the M-T2 variable, which combines the transverse mass of each lepton and the missing transverse momentum. No excess of events is found over the standard model predictions. Exclusion limits are placed at 95% confidence level on the production of top squarks up to masses of 208 GeV for models with a mass difference between the top squark and the lightest neutralino close to that of the top quark.Peer reviewe

Search for Higgs boson pair production in the gamma gamma b(b)over-bar final state in pp collisions at root s=13 TeV

A search is presented for the production of a pair of Higgs bosons, where one decays into two photons and the other one into a bottom quark-antiquark pair. The analysis is performed using proton-proton collision data at root s = 13 TeV recorded in 2016 by the CMS detector at the LHC, corresponding to an integrated luminosity of 35.9 fb(-1) . The results are in agreement with standard model (SM) predictions. In a search for resonant production, upper limits are set on the cross section for new spin-0 or spin-2 particles. For the SM-like nonresonant production hypothesis, the data exclude a product of cross section and branching fraction larger than 2.0 fb at 95% confidence level (CL), corresponding to about 24 times the SM prediction. Values of the effective Higgs boson self-coupling K X are constrained to be within the range -11 < K-lambda < 17 at 95% CL, assuming all other Higgs boson couplings are at their SM value. The constraints on K-lambda, are the most restrictive to date. (C) 2018 The Author(s). Published by Elsevier B.V.Peer reviewe

Search for an exotic decay of the Higgs boson to a pair of light pseudoscalars in the final state with two bquarks and two tau leptons in proton-proton collisions at root s=13 TeV The CMS Collaboration

A search for an exotic decay of the Higgs boson to a pair of light pseudoscalar bosons is performed for the first time in the final state with two b quarks and two tau leptons. The search is motivated in the context of models of physics beyond the standard model (SM), such as two Higgs doublet models extended with a complex scalar singlet (2HDM + S), which include the next-to-minimal supersymmetric SM (NMSSM). The results are based on a data set of proton-proton collisions corresponding to an integrated luminosity of 35.9 fb(-1), accumulated by the CMS experiment at the LHC in 2016 at a center-of-mass energy of 13 TeV. Masses of the pseudoscalar boson between 15 and 60 GeVare probed, and no excess of events above the SM expectation is observed. Upper limits between 3 and 12% are set on the branching fraction B(h -> aa -> 2 tau 2b) assuming the SM production of the Higgs boson. Upper limits are also set on the branching fraction of the Higgs boson to two light pseudoscalar bosons in different 2HDM + S scenarios. Assuming the SM production cross section for the Higgs boson, the upper limit on this quantity is as low as 20% for a mass of the pseudoscalar of 40 GeV in the NMSSM. (C) 2018 The Author(s). Published by Elsevier B.V.Peer reviewe

Measurement of inclusive very forward jet cross sections in proton-lead collisions at \sqrt{sNN} = 5:02 TeV

Measurements of differential cross sections for inclusive very forward jet production in proton-lead collisions as a function of jet energy are presented. The data were collected with the CMS experiment at the LHC in the laboratory pseudorapidity range −6.6 < η < −5.2. Asymmetric beam energies of 4 TeV for protons and 1.58 TeV per nucleon for Pb nuclei were used, corresponding to a center-of-mass energy per nucleon pair of \sqrt{sNN} = 5:02 TeV. Collisions with either the proton (p+Pb) or the ion (Pb+p) traveling towards the negative η hemisphere are studied. The jet cross sections are unfolded to stable-particle level cross sections with p_{T} ≳ 3 GeV, and compared to predictions from various Monte Carlo event generators. In addition, the cross section ratio of p+Pb and Pb+p data is presented. The results are discussed in terms of the saturation of gluon densities at low fractional parton momenta. None of the models under consideration describes all the data over the full jet-energy range and for all beam configurations. Discrepancies between the differential cross sections in data and model predictions of more than two orders of magnitude are observed

Search for the X(5568) State Decaying into B-s(0)pi(+/-) in Proton-Proton Collisions at root s=8 TeV

A search for resonancelike structures in the B-s(0)pi(+/-) invariant mass spectrum is performed using proton-proton collision data collected by the CMS experiment at the LHC at root s = 8 TeV, corresponding to an integrated luminosity of 19.7 fb(-1). The B-s(0) mesons are reconstructed in the decay chain B-s(0) -> J/Psi phi, with J/Psi -> mu(+) mu(-) and phi -> K+K-. The B-s(0)pi(+/-) invariant mass distribution shows no statistically significant peaks for different selection requirements on the reconstructed B-s(0) and pi(+/-) candidates. Upper limits are set on the relative production rates of the X(5568) and B-s(0) states times the branching fraction of the decay X(5568)(+/-) -> B-s(0)pi(+/-). In addition, upper limits are obtained as a function of the mass and the natural width of possible exotic states decaying into B-s(0)pi(+/-).Peer reviewe